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By performing quantum molecular dynamics (QMD) simulations, we investigate the equation of states, 
electrical and optical properties of the expanded beryllium at densities two to one-hundred lower than the 
normal solid density, and temperatures ranging from 5000 to 30000 K. With decreasing the density of Be, 
the optical response evolves from the one characteristic of a simple metal to the one of an atomic fluid. By 
fitting the optical conductivity spectra with the Drude-Smith model, it is found that the conducting 
electrons become localized at lower densities. In addition, the negative derivative of the electrical resistivity 
on temperature at density about eight lower than the normal solid density demonstrates that the metal to 
nonmetal transition takes place in the expanded Be. To interpret this transition, the electronic density of 
states is analyzed systematically. Furthermore, a direct comparison of the Rosseland opacity obtained by 
using QMD and the standard opacity code demonstrates that QMD provides a powerful tool to validate 
plasma models used in atomic physics approaches in the warm dense matter regime. 

Expanded liquid metals are focused in a series of important research areas because they evolve through 
various phases with associate transitions in the electrical and thermodynamic properties'. In this ther- 
modynamic regime, often labeled as warm dense matter (WDM), matter is partially dissociated, ionized and 
degenerate, which makes the description of the dynamical, electrical and optical properties extremely difficult. 
Experimental studies on electrical and thermodynamic properties of the metals with low critical temperatures 
such as mercury (Hg), rubidium (Rb) and cesium (Cs) have been widely made and the aspect of metal-nonmetal 
transition has been clarified"' *. The models to describe the metal-nonmetal transitions in very different physical 
systems from ordered solids to disordered fluids and plasmas and the important physical effects therein have 
already been discussed in detail by Mott, Hubbard, Anderson, Edwards et al. and Redmer et a/.'". It has been 
pointed out that the transitions of these very diverse systems are driven by the complex many-body nature of the 
respective interactions similarly". In the case of metals with higher critical points, the first experimental data were 
obtained on aluminum with the dynamic methods recently'''. Then, the metal-nonmetal transitions in expanded 
nickel'"*, titanium'*, iron'^, and tungsten"" were studied. One of the motivations of these studies was to find out the 
interrelation between the liquid- vapor and the metal-nonmetal transitions. However, to our knowledge, no data 
on expanded beryllium (Be), another refractory metal with high critical point, have been reported. As known, 
beryllium at warm dense regime is of particular importance because of its technological applications. Specifically, 
it is a potential ablator material in inertial confinement fusion (ICE) capsules and is also relevant for astrophysics. 
The optical and thermodynamic properties of Be decide the compressibility of the capsule, laser absorption, and 
instability growth at the fuel-ablator interface in ICE'^. 

The experimental and theoretical studies of warm dense Be mainly focused on its compressed state. Using the 
shock waves generated by underground nuclear explosives, Ragan'*" and Nellis et al. ''^ obtained the hugoniot of Be 
up to about 18 Mbar. The high-intensity laser in laboratory can even probe similar pressure range^". Recently, the 
dynamic structure factor in uncompressed and threefold compressed Be were calculated using the ab initio 
approach^'. Wang et al. reported the equation of state (EOS) and transport properties of Be for densities from 
4.0 to 6.0 g/cm' and temperatures from 1.0 to 10.0 eV based on the quantum molecular dynamics (QMD) 
simulations''^. However, the properties of thermal-expanded Be are still yet to be investigated. In the process 
of thermal expansion, a series of related questions show up: how the structure of Be evolves, whether the metal- 
nonmetal transition takes place, if yes, and then how does this transition affect its thermodynamic and optical 
properties. The answers to these questions would provide a good understanding of the structural, electronic and 



SCIENTIFIC REPORTS | 4:5898 | DDI: 1 0. 1 038/srep05898 



1 




10: T=6670+271*E 
12: T=6233+263*E 
18: T=5460+249*e" 
24: T=5605+235*E 
36: T=5377+222*E 



30 60 90 

Specific internal energy (kJ/g) 



120 



Figure 1 | (a) Pressure and (b) temperature versus specific internal energy 
along eight selected isochores ( V/Vq = 4 (square), 6 (circle), 8 (uptriangle), 
10 (downtriangle), 12 (diamond), 18 (star), 24 (sphere) and 36 (cross)). 
The lines are linear fits of the simulated results along the isochores 
correspondingly. 

dynamic properties of warm dense Be and the basic mechanisms 
that are responsible for these properties. From the theoretical 
point of view, QMD'^''^'' is particularly suited for this kind of pro- 
blems because no assumptions on the species involved or in the 
ionization state are needed. In particular, when combined with the 
Kubo-Greenwood formulation, QMD method can produce a con- 
sistent set of structural, electrical and optical properties from the 
same simulation. QMD has been successfully used on several 
expanded metals and predicted their metal-nonmetal transitions 
precisely'^""'^^^"^^. 

In this paper, we introduce QMD method, where the active elec- 
trons are treated in a fuU quantum mechanical way within the finite 
temperature density functional theory (FTDFT), to study the 
expanded Be. We calculate its EQS, electrical and optical properties 
to explore the evolution of Be as density and temperature change. The 
EOS are extracted from a series of QMD simulations of Be at densities 
of two to one-hundred lower than the normal solid density and 
temperature ranging from 5000 K to 30000 K. Using the 
configuration trajectories, the optical spectra are obtained from 
Kubo-Greenwood formula, from which the DC conductivity are 
determined. We analyze the variation of electronic density of states 
with respect to density and temperature to interpret the physics of the 
metal-nonmetal transition in expanded Be. The paper is organized as 
follows. In Sec. II, we briefly describe the simulation details and 
computational method in determining the optical conductivity. The 
results of EOS, electrical and optical properties are discussed in Sec. 
III. Finally, we close our paper with a summary of our main results. 

Results 

Equation of State. The EOS data for expanded Be are derived from 
the QMD simulations, which could describe various transient 
processes such as dissociation or association of chemical bonds 
and ionization or recombination induced by temperature or 
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Figure 2 | (a) Pressure and (b) specific internal energy versus the relative 
volume value along four selected isotherms. 

density without any approximation. The internal energy equals to 
the sum of the ion kinetic energy and the time average of DFT 
potential energy. The pressure consists of the contributions from 
the electronic and ionic components, which come from, 
respectively, the derivative of the total energy with respect to the 
Kohn-Sham electronic orbitals and the ideal gas expression since 
the ions move classically. In Fig. 1 the variations of pressure and 
temperature as a function of specific internal energy along eight 
selected isochores are shown. These relative volume values are 
y/Vo = 4, 6, 8, 10, 12, 18, 24, and 36, respectively, in which the 
normal solid volume Yq = 4.87 cmVmol. We can see from 
Fig. 1(a) that the EOS shows a systematic behavior for each 
isochore that the pressure and temperature increase almost linearly 
with the specific internal energy, similar as in expanded aluminum". 
The relations between temperature (in unit of K) and specific 
internal energy (in unit of kJ/g) are found to have the following 
forms of r = 9337 + 305£ for VIVq = 4, T = 8035 + 285£ for 
y/Vo = 6, r = 71 14 276£ for V/Vo = 8,T= 6670 + 271£ for V/Vo 
= 10, r = 6233 + 263£ for V/Vq = 12, T = 5460 + 249£ for V/Vo = 
18, r = 5605 + 235£ for V/Vo = 24, and T = 5377 + 222£ for V/Vq 
= 36. Since the experiment measured all quantities versus the 
internal energy variation while not temperature at present, it does 
not allow convenient comparisons with the other models. For 
expanded Be, the above relations map internal energies into 
temperatures, and could be used in future experimental study. In 
addition, we express the pressure (in unit of GPa) as a function of 
density p (in unit of g/cm"*) and temperature T as P = 4.83 — 
39.51454p -I- 0.00326pT through fitting all of these EOS data. It is 
indicated that the pressure is more dependent on density than 
temperature. Such a function should be applicable in hydrody- 
namic simulations for expanded Be conveniently. 

The density effects on the EOS are highlighted in Fig. 2, which 
presents the pressure and internal energy versus the relative volume 
value for four selected isotherms with T = 15000, 20000, 25000, and 
30000 K. As expected, the pressure decreases while specific internal 
energy increases with the relative volume systematically. In particu- 
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Figure 3 | The real part of frequency-dependent electrical conductivity for the four selected relative volumes of (a) VI Vq = 2; (b) VI Vq = 4; (c) VI Vq '■ 
8; and (d) V/Vo = 36 at temperatures of 10000 K (black), 20000 K (red) and 30000 K (blue). Data have been averaged over 10 to 15 uncorrelated 
QMD configurations. 



lar, it is noted that the internal energy keeps an almost linear relation 
with the relative volume until Be is expanded beyond eight times the 
normal solid volume. This means that the couplings between ions 
begin to become unimportant, indicating the appearance of an 
atomic phase. 

Electrical properties. In general, thermal-expansion of a metal could 
lead to metal- nonmetal transition. We examine the variation of 
electrical conductivity of Be over the extensive density- temperature 
ranges to clarify the metal-nonmetal transition. Using the Kubo- 
Greenwood formula, we calculate the real part of frequency- 
dependent electrical conductivity for points with relative volume 
ranging from 2 to 100 and temperature from 5000 to 30000 K. The 
real part of electrical conductivities o"i (to) for four selected relative 
volumes of VI Vq = 2, 4, 8, and 36 along the isotherms of 10000, 
20000, and 30000 K are presented in Fig. 3. It is obvious that these 
four sets of plots at different relative volumes exhibit different 
characteristics. In Fig. 3(a), the optical conductivity is in good 

agreement with a Drude form (o"£)(co) = — — -, with n* and t 

1 + or- 

respectively representing the conducting electron number density 
and relaxation time) with the maximum located at zero frequency, 
which indicates the nearly free-electron nature of the system. Such a 
behavior has also been reported for solid and three-compressed Be at 
high temperatures in Ref. 21. As expected, the dc conductivities 
become larger with increasing the density and temperature when 
comparing with our results, which predicts better conducting 
behavior. As the system is expanded to four times Vq, the optical 
conductivity loses the Drude character with the main peak shifted to 
about 1 eV, as displayed in Fig. 3(b). To describe this case more 
appropriately, the Drude-Smith modeP^ has been proposed in the 
following form 



0-05(01) = 



1 +03^Z^ 



1- 



c(l-«2T2y 



(1) 



which is a generalization of the classical Drude model through 
including the memory effects of the successive collisions. In the 
formula, c is a parameter of the memory effect. It equals to zero for 
the classical Drude model and falls into the range of 0 to — 1 when the 
main peak is shifted to higher frequency. The vanishing of dc 
conductivity could lead to the value of c = —1. With the relative 
volume increased to eight, the main peak continues to shift even to 
the energies larger than kgT, which corresponds to the case where the 
dc conductivity increases with temperature, opposite to the other two 
cases at lower relative volumes. It means that a metal-nonmetal 
transition takes place as Be is expanded to the volume larger than 
eight times Vg, which wiU be further discussed below. In addition, a 
small peak located between 3 ~ 4 eV seems to shows up, which 
becomes more evident at larger volumes. At the volume of thirty- 
six times Vo> an atomic transition line from 2s to 2p arises, with the 
location depending on the temperature. For example, the peak is at 
3.63 eV for T = 10000 K and 3.78 eV for T = 30000 K, in 
accordance with the atomic data. The shift comes from the 
increased ionization fraction of Be as the temperature is increased 
from 10000 K to 30000 K. It is confirmed that the expanded Be has 
evolved into an atomic phase from a condensed state. 

By fitting the optical conductivity spectra with the more 
general Drude-Smith model in the functional form of 



ffi(w) = 



(l + c)(l-hC02T2) 



C 1- 



(1-|-(«2t2) 



, we get the dc conduc- 



tivity (Tde> relaxation time t, and the parameter c. Meanwhile, the 
conducting electron density «* and effective ionization fraction Z 

can be estimated as rf = — - — and Z = , respectively, with 12 

being the volume of simulation cell and N, the number of ions in the 
cell. At the volume of two times of Vq, the fitting yields Z of about two 
and relaxation time of 1.57 X 10""" s for a temperature of 30000 K, 
in good agreement with an average ionization of 2.0 using both 
average atoms model""* and C0MPTRA04 program within the par- 
tially ionized plasma modeP^. This average ionization fraction leads 
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to a coupling parameter F = 14 and a degeneracy parameter of 0 = 
0.12, characteristic of a strong coupling and high degenerate plasma. 
At the largest volume of 100 times Vg and temperature of 20000 K, 
an analysis of the optical conductivity gives an ionization fraction of 
0.2, relaxation time of 5.37 fs and c = —0.725. The average atoms 
model'* gives an ionization fraction of 0.21, while COMPTRA04 
program"*^^ overestimates this value to be about 0.62. At such condi- 
tions, the system evolves into a weak coupled and low degenerated 
plasma with the corresponding coupling parameter and degeneracy 
parameter of 0.058 and 11.8, respectively. 

The Drude-Smith fitted conductivities are presented in Fig. 4. It 
illustrates two competing behaviors: a strong lowering of the elec- 
trical conductivity with decreasing density which may induce metal- 
nonmetal transition and a thermal activation of the conductivity due 
to ionization. The dc conductivity of Be along three isotherms of 
10000 K, 20000 K and 30000 K decreases with the density system- 
atically with a power-law behavior, and reaches a nearly constant 
value at the lowest density explored. Furthermore, we note that for 
volumes larger than eight times Vq, the enhancement of the dc con- 
ductivity with the temperature is clearly evident. Such a behavior 
could be seen clearer in Fig. 5. At the low relative volume values, 
the resistivity increases with the temperature, which is a character- 
istic of a metal. The sign of the temperature dependence of the 
resistivity becomes negative as Be is expanded to eight times Vq, 
demonstrating the behavior of a semiconductor. Such a metal-non- 
metal transition could also be confirmed from the point of the min- 
imum metallic conductivity of about 2000 fl^'cm"' for disordered 
systems"", because we can see that the dc conductivity of Be is smaller 
than this value before the crossing points of these three plots. 

A typical way to characterize the nonmetal-metal transition in 
expanded Be is the variation of electronic density of states (EDOS). 
In Fig. 6 we give the averaged EDOS at 10000 K for four different 
densities with a broadening of 0.03 eV. It can be seen that at V/Vg = 
2, EDOS clearly exhibits a free-electron-like behavior. As the volume 
increases, the EDOS is depressed in the vicinity of the Fermi level. At 
the largest volume of thirty-six times Vg presented in Fig. 6(d), the 
EDOS shows a local atomic behavior and a gap opens. 

Optical properties. The evolution of the media with lowering the 
density becomes clearer by examining the variation of the optical 
absorption coefficient along the 10000 K isotherm, as shown in 
Fig. 7. It is observed that the absorption coefficient increases 
rapidly in the photon energy range of a few eV, which mainly 
comes from the free-free contribution (also called inverse 
bremsstrahlung). As the volume is expanded, the most prominent 
change in the absorption coefficient is the shape and location of its 
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Figure 4 | The variation of the dc conductivity of Be as a function of 
density for three isotherms. 
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Figure 5 | The resistivity variation with respect to the temperature along 
the five selected isochores of V/ Vq = 1 (square), V/Vq = 2 (circle), V/ Vg = 
4 (uptriangle), V/Vq = 8 (downtriangle), and V/Vq = 12 (diamond). 

main peak. At the lowest volume, the absorption coefficient shows 
the characteristic of a simple metal as described by the Drude model, 
with the maximum corresponding to a collective behavior associated 
with the transitions to the lowest excited states. With the volume 
further expanded up to twelve times Vq, more local features arise, 
especially the atomic transition line. At the volume of 100 times Vg, 
the peak become sharper due to the decrease of the ionization 
fraction of Be. Its location at about 3.63 eV corresponding to the 
position of 2s to 2p atomic transition. In addition, a valley around 
2.0 eV is associated with the opening of the band gap as discussed 
above. 

In Fig. 8, we inspect the variation of the index of refraction with 
photon energy for different densities. For lower volumes than eight 
time Vq, the index of refraction decreases smoothly down to a min- 
imum at an energy corresponding to the plasma frequency. As the 
volume increased to twelve time Vq, a mild maximum at a photon 
energy about 3.58 eV arises, followed by a marked minimum shifted 
to lower energy compared to that for lower volume. This suggests 
that atomic phase begins to appear under such density-temperature 
conditions. With the volume further expanded, the index of refrac- 
tion first decreases to a minimum corresponding to a plasma fre- 
quency and then exhibits resonant profile at the position of atomic 
transition line, which is featured uniquely by the noninteracting 
atoms or ions. In fact, the plasma frequency could also be evaluated 

using the formula co = with the known ionization fraction. 

me 

For example, at the lowest volume of two times Vq, the plasma 
frequency is calculated to be 1 3 .06 eV corresponding to an ionization 
fraction Z = 2, while for the highest volume of 100 times Vq, the 
plasma frequency is obtained as 0.21 with an ionization fraction of 
0.027, in agreement with the behavior of the index of refraction 
shown in Fig. 8. Therefore, QMD is proved useful to describe the 
variation of the optical properties for a wide range of conditions from 
isolated atomic to the warm dense matter regime. 

Discussion 

To understand the variations of the properties of expanded Be, we 
inspect the differences in the Rosseland mean opacities obtained 
using QMD and the standard opacity code LEDCOP". The 
Rosseland mean opacity is calculated using Eq. (2) and based on 
QMD simulations for different density-temperature conditions. 
Figure 9 displays the variation of the Rosseland mean opacity for 
Be as obtained for a fixed temperature (10000 K) with increasing 
density (left panel) and fixed densities {V/Vg = 4 and V/Vq =12) 
but increasing temperature (right panel). Comparisons with the data 
from LEDCOP are also included. Different from the QMD method 
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Figure 6 | Electronic density of states of expanded Be at 10000 K for four densities. Data have been averaged over 10 to 15 uncorrelated configurations. 
The zero of the energy scale shows the position of the Fermi level. 



which gives a global optical spectra, LEDCOP opacity code calculates 
the contributions from atomic photoionization (bound-free), inverse 
bremsstrahlung (free-free), scattering and bound-bound transition 
separately. As seen in Fig. 9(a), for low densities, LEDCOP and the 
QMD show good agreement for Rosseland mean opacity. However, 
as the density increases, the LEDCOP and the QMD calculations 
depart, which results from the fact that various density effects, such 
as pressure ionization, are only introduced in a phenomenological 
fashion in LEDCOP. In addition, it is noted that the trend of the 
variation of QMD Rosseland mean opacity with expanding the 
volume changes beyond eight times Vq, which just coincides with 
the expansion induced metal-nonmetal transition. In Fig. 9(b), the 
Rosseland mean opacity decreases with temperature for these two 
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Figure 7 | Absorption coefficient a (to) as a function of photon energy at 
T = 10000 K for different densities. Bar at 4fcBT represents maximum 
region of contribution to the Rosseland mean opacity. 



densities, because at temperature of a few eV the opacity mainly 
comes from the free-free contribution. 

In summary, we have systematically investigated the equation of 
state, electrical and optical properties of expanded Be at densities two 
to one-hundred lower than the normal solid density and tempera- 
tures ranging from 5000 to 30000 K through performing QMD 
simulations. The optical spectra demonstrate that as the density 
decreases, the system evolves from a simple metal to an atomic fluid. 
Through analyzing the variation of electrical conductivity, it is found 
that metal-nonmetal transition takes place in expanded Be at density 
about eight lower than the normal solid density, with the conducting 
electrons becoming localized. Furthermore, it is proved that QMD 
provides a powerful tool to describe from warm dense matter to 
atomic fluid. 
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Figure 8 | The index of refraction u (») as a function of photon energy at 
r = 10000 K for different densities. 



SCIENTIFIC REPORTS | 4:5898 | DDI: 1 0. 1 038/srep05898 



5 



1.5x10 
l.Oxlo' 



5.0x10 



0.0 




QMD T=10000K 
LEDCOP T=11600K 



10 





(b) _ 


-□- QMD V/Vo=4^ 




-O-QMD V/Vo=12 




10000 20000 


30000 



1.5x10 



l.OxlO' ^ 



0.0 



V/Vo 



Temperature (K) 



Figure 9 Variation of QMD and LEDCOP Rosseland mean opacities as a function of relative volume for a fixed temperature (left panel), and as a 
function of temperature for two fixed densities (right panel). 



Methods 

In the present application, ionic structures are first generated with the QMD package 
VASP {Vienna ah initio simulation package) developed at the Technical University of 
Vienna^^'^^, in the framework of FT -DFT^ The interactions between ions and 
valence electrons are described using the all-electron projector augmented wave 
(PAW) potentials*^ The exchange- correlation energy is treated on the level of 
generalized gradient approximation (GGA) using the Perdew-Wang 91 parametri- 
zation'*\ The electronic states are populated through a Fermi-Dirac distribution, with 
electronic temperature set equal to that of ions. Ions move classically according to the 
forces, which originate from the electron density and the ion-ion repulsion. 

All of our QMD simulations employ the canonical (NVT) ensemble with Nose- 
Hoover thermostat"*^. We select the densities ranging from two to one-hundred lower 
than the normal solid density (V/Vq — 2 — 100) and temperatures from 5000 to 
30000 K which cover a broad phases. Simulations are done at F point with 54 atoms 
for the densities from two to eighteen lower than the normal solid density and 8 atoms 
for the densities from twenty-four to one-hundred lower than the normal solid 
density in a cubic cell of the size appropriate to the considered density. More atoms 
are tested but the difference could be negligible. The plane-wave cutoff is fixed at 
500.0 eV, which is tested to give good convergence for both total energy and pressure. 
For electronic properties calculations, 4X4X4 Monkhorst-Pack^'^ scheme k points 
are used to sample the Brillouin zone. Integration of the equations of motion proceed 
with time step of 0.5 — 1.0 fs for different pressure- temperature ranges. After about 
3 ps the system is equilibrated and the subsequent 5 ps are taken to calculate the EOS 
and electronic properties as running averages. The ion temperature T, is fixed using 
Nose-Hoover thermostat"*^. Additional convergence tests with respect to particle 
number and k-point Brillouin sampling are performed for several pressure-temper- 
ature conditions. In all cases, the internal energy is converged to less than 0.01 eV/ 
atom and the pressure within an error bar of 2%. 

Based on QMD simulated trajectories, 10 — 15 configurations spaced by at least the 
correlation time are selected. For each of these configurations, the optical conduct- 
ivities are obtained with the Kubo- Greenwood formulation^^*^ {atomic units), 

9„ N N 3 

k ;=1 i=l a = l (2) 

xK'i';.k|V,|1',k>|',5(9.k-c,,k-fim), 

where oj is the frequency, Q the atomic volume, and Wthe total number of bands used. 

and f ,.k are the electronic eigenstates and eigenvalues for the electronic state i at a 
given k point with the w {k) k — point weighting factor, /(c/k) stands for the Fermi 
distribution function. Other optical coefficients can be directly derived from the 
optical conductivity. In addition, the static limit w ^ 0 of optical conductivity gives 
the dc conductivity a^^. As a particular application, the Rosseland mean opacity^^ is 
commonly used to represent the radiative property, which is expressed as 

"r(P^'^) dcon^co)B{mJ) ' ^ ' 

where n (co) and y. {co) are the frequency- dependent index of refraction and 
absorption coefficient, respectively. B' (co, T) represents the derivative of the nor- 
malized Planck distribution with respect to temperature at a given T, a slowly varying 
function which peaks around 4fcaT. 
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